Vibrational properties of inclusion complexes: the case of indomethacin-cyclodextrin 
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Vibrational properties of inclusion complexes with cyclodextrins are studied by means of Raman 
spectroscopy and numerical simulation. In particular, Raman spectra of the non-steroidal, anti- 
inflammatory drug indomethacin undergo notable changes in the energy range between 1600 and 
1700 cm -1 when inclusion complexes with cyclodextrins are formed. By using both ab initio quan- 
tum chemical calculations and molecular dynamics, we studied how to relate such changes to the 
geometry of the inclusion process, disentangling single-molecule effects, from changes in the solid 
state structure or dimerization processes. 

PACS numbers: 61.43.-j, 61.43.Fs, 63.50.+X 
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I. INTRODUCTION 

Inclusion complexes with cyclodextrins have been lit- 
tle studied by means of Raman spectroscopy, although 
the usefulness of this technique in such systems is now 
acknowledged^^ Cyclodextrins (CD) are a family of nat- 
ural or synthetically modified cyclic molecules, consisting 
typically of six (cv-CD), seven (/3-CD) or eight (7-CD) 
glucopyranose units. In water, they take on the pecu- 
liar 3-dimensional structure of a truncated cone, with a 
slightly soluble outer surface and a hydrophobic central 
cavity, whose size depends on the number of glucose units 
of the molecule. A remarkable property of CD in aque- 
ous solution, is their ability to form host-guest inclusion 
complexes with a wide variety of organic and inorganic 
molecules, provided that the guest molecule is less po- 
lar than water* 3 In particular, inclusion complexes of CD 
with non-polar drugs are a topic of current interest, be- 
cause these non-covalent complexes increase not only the 
aqueous solubility of drugs, but also their chemical sta- 
bility and bioavailability. 

Such complexes have been usually investigated in 
aqueous solution by ultra-violet and visible absorption 
spectroscopy, fluorescence techniques, and NMR spec- 
troscopy, while X-ray and neutron diffraction are typi- 
cally utilized in the solid state, but the latter techniques 
(unlike Raman) require crystalline samples. The effect 
of the inclusion process on the guest molecules has also 
been the subject of some Raman studies^*£*& in particu- 
lar, Raman and FT-IR investigations of the interactions 
between CD and some non-steroidal anti-inflammatory 
agents (like diclofenac sodium^ ketoprofen£ and amor- 
phous piroxicam 9 ), have shown significant changes in the 
vibrational spectrum of the complexed guest molecules 



with respect to the free ones. 

Among the non-steroidal anti-inflammatory drugs, in- 
domethacin (IMC, Fig. (IJ is widely used as an analgesic 
drug in the treatment of rheumatoid arthritis, as well as 
in other degenerative joint diseases, and recently it has 
also shown anti-tumor activity, with an important role 
in the inhibition of human-colon adenocarcinoma cell- 
proliferationiiii The molecule of IMC has been charac- 
terized by different spectroscopic techniques, and it has 
been observed that, in the crystalline state, it exhibits 
three polymorphic forms (called a, (3 and 7, respectively), 
depending on the sample history. 11 Due to its chemical 
structure, IMC is poorly soluble in water, resulting in 
high local drug concentrations responsible for ulcerations 
in the tissues, and this reduces its therapeutic applica- 
tions. A strategy to affect the solubility and chemical 
stability of IMC in water, which is actually employed in 
commercial drugs, consists in the preparation of inclusion 
complexes with CD. 

Although some models were proposed for the arrange- 
ment of the guest molecule into the cavit y 12 ! 13 unequivo- 
cal conclusions were not reached. Moreover, although the 
vibrational spectra of the various forms of free IMC have 
been widely studied in the literature by FT-IR and Ra- 
man spectroscopy^ no Raman data have been reported 
on its inclusion complex with CD. 

In the present work we discuss the results of Raman 
scattering experiments and numeric simulations on the 
IMC-CD inclusion complexes in solid state, which pro- 
vide new insight into the structure of the complexes, and 
into the effect of the inclusion process on the guest. Ad- 
ditional information is obtained by the comparison of 
these results with those relative to (i) 5-metoxyindol-3- 
acetic acid, (ii) the inclusion complexes formed by CD 
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FIG. 1: Chemical structure of indomethacin, with numbering 
of atoms adopted for the molecular dynamics computations. 



with 4-cloi'obenzoic acid and (iii) with the sodium salt of 
IMC (NalMC). Additional experiments of electrospray- 
ionization mass spectrometry and NMR were also per- 
formed. 



II. EXPERIMENTAL AND COMPUTATIONAL 
DETAILS 

All chemical compounds were purchased from Sigma- 
Aldrich and used without further purification. Raman 
spectra were recorded using a Jobin-Yvon HR800 micro- 
Raman confocal system with a 100 x objective, equipped 
with a liquid-nitrogen-cooled CCD detector. Exciting ra- 
diation at 632.8 nm (20 mW) was provided by a He-Ne 
laser. The Raman spectra were generally recorded using 
dried samples. In some cases, when the amount of avail- 
able complex was small, in order to improve the Raman 
signal, Surface Enhanced Raman Scattering (SERS) was 
used. In these cases, drops of liquid sample were repeat- 
edly deposited on a suitable metallic substrate (Ti-6A1- 
4V alloy) and dried by evaporation, thereby increasing 
the complex concentrationii£ii£ii£ Mass spectrometry ex- 
periments were performed on a Bruker Esquire-LC™ ion 
trap mass spectrometer equipped with an atmospheric- 
pressure Electrospray ionization (ESI) source, in the volt- 
age range 99-180 V; the total ion current was acquired 
with a scan range m/z 100-2000. The sample was di- 
rectly infused as a 1 : 1 MeOH/H^O solution at a flow 
rate of 2 /xl/min into the ESI source in either positive- or 
negative-ion mode. 1 HNMR spectra were taken at room 
temperature with an Avance 400 Bruker spectrometer; 
1 H was recorded at 400 MHz, with S values in ppm in 
D 2 relative to the HDO residual signal (5h 4.70 ppm) 
and J values in Hz. 

All the ab initio quantum chemical calculations were 
performed with the GAUSSIAN 03 program suited uti- 



lizing unrestricted DFT 19 and the non-local B3LYP func- 
tional hybrid method was employed^ as well as un- 
restricted Hartree-Fock (UHF). The standard 6-31G(d) 
basis se1*2I was used for the geometry optimization and 
frequency analysis of free indomethacin base pairs. Vi- 
brational frequencies were computed at the same level of 
theory, as well as IR and Raman spectra. In addition, 
we have also obtained the amplitude of the displacements 
of the atoms for the different normal modes. When the 
vibrational frequencies have been computed, the scaling 
factor of 0.962S was taking into account; this has to be 
introduced to correct for anharmonicity, as well as for 
the systematic errors inherently present in the calculated 
harmonic contributions^. 

Molecular dynamics simulations were performed 
by utilizing the Amber packaged* In particular, the 
molecular dynamics runs have been performed via the 
module sander, while the energy minimization and the 
search for eigenvalues and eigenvectors of the dynamical 
matrix have been obtained via the module nmode. 
The force field employed is GAFF (Generalized Amber 
force field3£,). In some runs, the solvent effects have 
been incorporated explicitly, by including TIP3 water 
molecules^ The starting structures of IMC and /3CD 
have been obtained from the Protein Data Banki2£ 
The coordinates for hydrogen atoms were generated by 
means of the WWW PRODRG server^ Partial charges 
were computed by means of the module antechamber of 
the Amber package, utilizing the AM1-BCC model^ 

The IMC-CD complex was prepared according to the 
procedure adapted from the literature^ /3CD (50.3 mg, 
44.4 ^mol) was dissolved in water (4.2 ml) obtaining 
a 10.6 /xM solution, which was added to an equimolar 
amount of dry IMC (15.9 mg, 44.4 /xmol). The resulting 
dispersion was stirred at 60° C for 2 hours and cooled at 
room temperature. After 3 hours, the presence of the 
complex was indicated by the precipitation of a white 
solid, which was collected by drawing the liquid phase 
with a syringe, dried in a vacuum chamber (over P2O5 
as dehydratant), and finally powdered. The complexes 
between hydroxypropyl /3CD (HP/3CD) and IMC, and 
between /3CD and 4-clorobenzoic acid or sodium salt of 
IMC, were obtained under the same procedure. 

By 1 HNMR analysis in D2O on the complexes of both 
IMC and NalMC with /3CD or HP/3CD, we have ob- 
tained the same *H chemical shifts previously reported 
by Mylesi^ and Fronza^ 



III. EXPERIMENTAL RESULTS 

Before the Raman scattering measurements, prelimi- 
nary experiments were performed on the complexes using 
mass spectrometry. In fact, when applied to the analy- 
sis of CD complexes, the negative ion ESI-MS technique 
provides information on the stoichiometry of the com- 
plex in gas-phase. In the case of the IMC-/3CD complex, 
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FIG. 2: ESI(-)-MS spectrum for the IMC-/3CD complex. 



as shown in Fig. [3 the cluster signals of the pseudo- 



molecular ion [M comp i ex -H] ~ of the 1:1 complex is de- 
tected at m/z— 1490, and the signal of the ion [M<7d-H]~ 
of free /3CD at m/z— 1133. No signals arising from 1 : 2 
complex were observed. Tandem MS/MS experiments, 
carried out in the ion trap through isolation of the clus- 
ter of the 1 : 1 complex at 1490 m/z, and CID fragmen- 
tation analysis, yielded the daughter ion at m/z 1133 
by neutral loss of IMC. Similar behavior was observed 
for IMC-HP/3CD and IMC sodium salt-/3CD complexes, 
supporting 1:1 stoichiometry of the complexes in these 
cases as well. 

In Fig. |21 we report the Raman spectra of solid 7-IMC 
(a), 4-chlorobenzoic acid (b), 5-metoxyindol-3- acetic acid 
(c), and the 1:10 physical mixture of IMC and /3CD (d), 
in the energy range 1500-1750 cm" 1 . The two latter 
molecules are of interest in the present study because 
their chemical structures are very similar to subunits of 
IMC, as shown in Fig.0] 
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FIG. 3: Top: Raman spectra of indomethacin (a), 4- 
chlorobenzoic acid (b), 5-metoxyindol-3-acetic acid (c), and 
1:10 physical mixture of IMC and /3CD (d). Throughout the 
paper, different spectra are vertically shifted for the sake of 
clarity. 
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FIG. 4: Chemical structure of 5-metoxyindol-3-acetic acid 
(1), and 4-chlorobenzoic acid (2). 



We focused on the this energy range because it contains 
no free-/3CD peaks, and complexation-induced changes 
on IMC spectrum can be readily observed. Moreover, 
this is also the range where C=0 stretching energies are 
expected, so that combining experiment and simulations, 
useful information on the geometry of the complex (for 
instance, the position of the involved CO group) may be 
obtained from the features of Raman spectrum. IMC 
contains both a C=0 amide group (atoms N28-C29-O30 
in Fig. nj and a C=0 carboxylic group (atoms C17- 
O18-O19-H20), which are spatially well separated in the 
molecule. The comparison among the Raman spectra of 
Fig. |3| seems to indicate that, actually, the peak of IMC 
at 1698 cm -1 corresponds to the amide C=0 stretch. In 
fact, this peak does not appear in the spectra of com- 
pounds (1) and (2) of Fig. QJ, which do not contain the 
amide group; it is also be noted, that this assignment is in 
agreement with Taylor et al. 4^ Moreover, in the following 
we shall see that ab initio quantum chemical calculations 
and molecular dynamics computations confirm this con- 
clusion. Thus we may focus on the peak at 1698 cm -1 
(from now on the corresponding frequency will be la- 
belled lucoNi while the frequency of the carboxylic group 
will be labelled locooh), and on its changes due to the in- 
clusion process, which are shown in Fig. [51 By comparing 
the spectrum of free IMC (a) to those of the complexes 
formed by HP/3CD (b) and /3CD (c), we note a marked 
broadening of the peak at ujco n , as well as a shift from ss 
1700 cm -1 to « 1670cm -1 . Such a large shift for lucon 
(~ 30 cm -1 ) is rather unusual in the study of the com- 
plexation processes where, generally, smaller variations 
are observedi&S From the above results one might infer 
that it is the amide C=0 group of IMC (and the neigh- 
boring atoms) to be mainly affected by the inclusion pro- 
cess. However, in order to understand the precise way in 
which such atoms are influenced by complexation, other 
possible effects should be discussed. Generally speak- 
ing, changes of vibrational frequencies might be due to 
the presence of the uncomplexed guest in dimeric form. 
This is, for example, what happens to the 4-chlorobenzoic 
acid which, in the crystalline state, may exist also in 
dimeric form due to intermolecular hydrogen-bonded in- 
teractions, and whose Raman spectrum is shown in Fig. [21 
(b): here, the broad peak at 1628 cm -1 has been assigned 
to the carboxylic C=0 stretch of the dimerj^i and thus 
corresponds to ujcooh- In Fig. we report the spectra of 
free 4-chlorobenzoic acid (a), of the 1 : 10 physical mix- 
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FIG. 5: Raman spectra of free IMC (a), IMC-HP/3CD com- 
plex (b), IMC-/3CD complex (c). 

ture of 4-chlorobenzoic acid-/3CD, respectively (b), and 
of the /3CD-complexed 4-chlorobenzoic acid (c). In case 
(b), only 4-chlorobenzoic acid monomers are expected: 
as a matter of fact, the breakdown of hydrogen bond- 
ing patterns of 4-chlorobenzoic dimers is observed in the 
shift of the C=0 stretch from 1628 cm" 1 (Fig.^a)) to 
higher energy (w 1684 cm' 1 ) (Fig.HJb)). Moreover, the 
comparison between Fig. (b) and (c) shows that the 
complexation process further shifts the peak correspond- 
ing to C=0 stretch by ~ 10 cm -1 . Note that this high- 
energy shift for the 4-chlorobcnzoic acid is at odds with 
the IMC case, where a decrease of locon was observed 
upon complexation. 
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FIG. 6: Raman spectra of free 4-chlorobenzoic acid (a); 1:10 
physical mixture of 4-chlorobenzoic acid-/3CD, respectively 
(b); 4-chlorobenzoic acid-/3CD inclusion complex (c). 

The case of the 4-chlorobenzoic acid shows that com- 
plexation may affect the Raman peaks in at least two 
ways: (i) by breaking the intermolecular bonds, (ii) by 
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FIG. 7: Raman spectra of free IMC sodium salt (a), and IMC 
sodium salt-/3CD complex (b). 



modifying the strength of the interatomic bonds (C=0 
in this case). In order to check the possibility that the 
former effect may affect the complexed form of IMC, we 
have measured the Raman spectrum of a 1 : 10 phys- 
ical mixture of IMC-/3CD (Fig. 03(d)). In the energy 
range between 1500 and 1750 cm -1 , this spectrum ex- 
hibits the same characteristic peaks as in Fig. (a), indi- 
cating that the shift of the amide C=0 stretch, observed 
in Fig. [5] (a) and (b) , is actually related to complexation 
and not to the cleavage of hydrogen bonding patterns of 
IMC dimers. 

Similar results have been obtained on the inclusion 
complex formed by /3CD with IMC sodium salt (Fig. [TJ: 
also in this case, we clearly observe a broadening and shift 
to lower energy of the peak at 1675 cm -1 , which is gen- 
erally assigned to the amide C=0 stretching, when the 
guest is included in /?CD. All these results confirm the 
trend already observed in the vibrational spectra of IMC 
complexes, and further support the hypothesis that the 
amide C=0 of IMC is directly involved in the process of 
complexation with /3CD. Finally, the vibrational prop- 
erties of IMC-/3CD inclusion complexes have been also 
investigated in the low-frequency region (10-100 cm -1 ), 
using a Jobin Yvon U1000 Raman spectrometer with Ar- 
laser excitation at 514.5 nm; preliminary results seem to 
evidence, in this spectral region as well, the presence of 
significant differences, probably due to extended vibra- 
tions of the system^ 

Actually, there is also a third phenomenon that might 
induce changes in the experimental spectra, namely the 
possible change of the solid state of IMC due to the com- 
plexation. In general, below its glass transition temper- 
ature IMC crystallizes in the stable 7-form, whereas at 
higher temperatures the formation of the metastable a- 
form is predominant »Ai The polymorphic a and 7 forms 
differ not only in the packing density of the crystals, but 
also in the unit cell arrangements^ leading to a deep 



difference in their Raman spectral Thus, one could ar- 
gue that the large structural change undergone by the 
solid state due to complexation, could be related to the 
observed shift of ujcon- In order to investigate this pos- 
sibility, we performed numeric simulations describing the 
complexation process between a single molecule of IMC 
and a single molecule of (3CD. We shall see that they 
reproduce sufficiently well the experimental findings and 
that structural changes do not need to be taken into ac- 
count. 




IV. NUMERIC SIMULATIONS 

A. Ab initio quantum chemical computation 

Focusing only on the energy range described in the 
experimental section, the main result from the ab initio 
computation is that, in the region 1500 — 1750 cm -1 (a 
complete assignment of frequencies for free IMC has been 
attempted in Ref. 3 i), the eigenvector with the largest 
Raman signal involves the largest displacement over the 
atoms C29-O30, i.e. it corresponds to the amide C=0 
stretching. The value for u>con found from the UHF ab 
initio computation (1696 cm' 1 ) is rather close to the ex- 
perimental one, as shown in Fig. |SJ Furthermore, in the 
comparison between experimental and computed Raman 
intensities, the signals falling in the region between about 
1600 and 1630 cm -1 can probably be assigned to aro- 
matic C=C stretch bond of IMC. The vibrational analy- 
sis of complexes has been performed using classical sim- 
ulation (see next section) , which allows larger systems to 
be simulated with respect to the ab initio calculation. 
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FIG. 8: Full line: experimental Raman spectrum of free IMC; 
Full circles: Raman intensities as computed by UHF. 



FIG. 9: Starting configuration of the IMC-/3CD complex uti- 
lized in cases (3) and (4) (see text). In the former case, 
molecules of water have been added. 



B. Molecular dynamics computation 

The advantage of numeric simulations, with respect to 
the experiments, is that the positions of the IMC and 
CD atoms can be located exactly. On the other hand, 
numeric simulations of the complex are rather CPU-time 
demanding, especially when the effects of the solvent are 
taken into account. The length of our runs ranged from 
1 ms, for a molecule of IMC in vacuum, to only a few 
ns when thousands of molecules of water were added. To 
answer in a reliable way the question about the geometry 
of the inclusion process, one should perform long enough 
runs of many IMC and /3CD molecules, in the presence of 
a solvent, in order to detect many spontaneous inclusion 
events, and then perform a statistical study. As we shall 
see, because of the long CPU-time needed, this goal has 
been only marginally achieved, since only one of such 
events has been observed. On the other hand, it has been 
possible to gain further physical insight by supporting 
this result with somewhat less ambitious simulations. 

In order to study the inclusion of IMC in /3CD we 
have simulated the following systems: (1) single IMC 
molecule in vacuum, (2) single IMC molecule in a box of 
1114 molecules of water, (3) single IMC-/3CD inclusion 
complex in vacuum, (4) single inclusion complex in a box 
of 633 molecules of water and (5) 9 IMC molecules plus 
9 /3CD molecules in a box of 1400 molecules of water. 
In systems (3) and (4) the inclusion complex was pre- 
pared by docking the CI atom of IMC into the cavity of 
/3CD, and minimizing the energy of the resulting configu- 
ration (see Fig.|5J|. During the molecular dynamics runs, 
we recorded the distance d between the CI atom of IMC 
and the center of mass of /3CD, in order to determine 
the formation of an inclusion complex. More precisely, a 
complex is defined by the fact that d ~ (or, at least, 
much lesser than the typical distance between free IMC 
and /3CD). Moreover, starting from the thermodynamic 
equilibrium configuration, we sought for the minima of 
the energy potential. In the minimum configurations we 
computed the eigenvalues u> and the corresponding eigen- 
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vectors e,-(u;) of the dynamical matrix, in order to deter- 
mine the vibrational properties of the different systems. 
From such quantities we have obtained the density of vi- 
brational states g(ui) and the quantity Yi = e*j • efj, which 
measures the amplitude of the vibration of the i-th atom 
due to the mode of frequency uj. The thermodynamic 
parameters chosen for the simulations are room temper- 
ature and room pressure, with the coupling to the ex- 
ternal bath introduced via a weak coupling schemed In 



system (5) 
system (4) 
system (3) 
N.N. IMC-CD distance 




100 
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FIG. 10: Distance d between the CI atom of IMC, and the 
center of mass of /3CD, as a function of time for different 
simulated systems (see text). 

Fig. (|10fl we show the behavior of the distance d during 
different molecular dynamics runs. The solid line repre- 
sents the distance between the CI atom of IMC and the 
center of mass of /3CD, for a given pair chosen among the 
9 IMC plus 9 /3CD molecules of system (5). After 1 ns, 
d reaches a value of approximately 2 A, which we inter- 
pret as the signature of an inclusion event, with the CI 
located quite close to the center of mass of /3CD. In our 
runs, we did not observe other inclusion processes, thus 
we are not able to perform a detailed study of the phe- 
nomenon (characteristic times, etc.). It is interesting 
to note that d is much smaller than the position of the 
first peak of the pair correlation function between free 
IMC and CD molecules (horizontal line). Furthermore, 
the value of d reached in the spontaneous inclusion pro- 
cess found in system (5) is very similar to the stationary 
value of the distance in the run of system (4) (dashed 
line), where the complex was initially created and put in 
a box of water. Both results support our claim that in 
system (5) (which is the most similar to an experimen- 
tal sample) we actually observed a single complexation 
process. Interestingly, the complex is not stable without 
water: in fact, we see that although system (3) (dotted 
line) starts from the same initial configuration as system 
(4), the absence of water leads to an immediate destruc- 
tion of the complex. 

Unlike the ab initio quantum chemical computations 
discussed previously, from classical molecular dynamics 
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FIG. 11: Top: density of vibrational states g(u>) for the sim- 
ulated free IMC (solid line) and /3CD (dashed line). Bot- 
tom: projection of the eigenvectors corresponding to ujcon 
and uicooh over the atoms of IMC, labelled as in Fig. 1. 
Note that in both cases there are also large contributions from 
some hydrogen atom. 



(MD) simulations it is not possible to obtain in a straight- 
forward way the Raman spectra; in fact, the Raman 
signal I(lu) is cx g(u) C(cu), where C(ui) is the so-called 
coupling function which, in classical computations, is un- 
known unless the mechanism of polarizability modulation 
is specified^ On the other hand, with MD simulations 
one can obtain the eigenfrequencies and the eigenvectors 
of the inclusion complex. 

The interpretation of the density of vibrational states 
gibj) proceeds through different steps. In Fig.^]we iden- 
tify the frequencies ujcon— 1628 cm -1 (Raman active ac- 
cording to the ab initio calculation, see Section IV A) and 
ucooh= 1649 cm -1 (Raman inactive), as the ones cor- 
responding to the modes which are most localized on the 
atoms of the amide and of the carboxylic C=0 groups, 
respectively. There is a slight difference between the val- 
ues of these two characteristic frequencies and the ones 
computed with the ab initio study, and between ujcon a s 
computed with the MD approach and the experimental 
value. These discrepancies might be ascribed to the force 
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field utilized. 
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FIG. 12: Top: projection of the eigenvector corresponding to 
ucon over the atoms of IMC (labelled as in Fig. 1), in free 
IMC (solid line, system (1)) and IMC-/3CD complex (dotted 
line, system (4)). Bottom: projection of the eigenvector cor- 
responding to ujcooh in free IMC (solid line, system (1)) and 
IMC-/3CD complex (dotted line, system (4)). 

In Fig. E| is shown the effect of complexation on the 
eigenvectors corresponding to locon and ujcooh- In 
fact, in the complexed system (4), we have looked for the 
frequencies whose corresponding eigenvectors were the 
most similar to the ones of free IMC shown in Fig. ^2 
While ujcooh is not affected by the inclusion process, 
one can see that ujcon, corresponding to the largest ex- 
perimental Raman signal, shifts to lower frequencies by 
~ 14 cm -1 . Thus, from a qualitative point of view, we 
find results in agreement with the experimental findings 
described in the above section. Quantitatively, one has 
to note that this shift is about one half of the experimen- 
tal one. It seems reasonable to ascribe this difference to 
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FIG. 13: Density of vibrational states g(uj) for free /3CD 
(solid line), for 4-chlorobenzoic acid (dashed line), and for 
4-chlorobenzoic-/3CD-complex (dotted line). 

the limitations of the GAFF force field in describing a 
molecule of the complexity of IMC, rather than to the 
fact that in these simulations the structural properties 
of the solid state (and their possible modifications due 
to the complexation) are not taken into account. As a 
matter of fact, in Fig. 1131 we show that for a much sim- 
pler molecule like the 4-chlorobenzoic acid (for which the 
GAFF force field is expected to be more reliable), the 
MD simulations show a shift of ujcooh to higher en- 
ergy of ~ 10 cm -1 , which is in good agreement with 
the one found in experiments in the complexation of the 
monomcric form of 4-chlorobenzoic acid. 



V. CONCLUSIONS 

In the present paper we have reported the results of 
a multi-technique study of the complexes formed by CD 
with the molecule of IMC. In particular, we have per- 
formed Raman scattering experiments and we have com- 
pared the experimental findings with the outcome of ab 
initio quantum chemical calculations and of molecular 
dynamics simulations. Focusing on the energy range 
1500 — 1700 cm -1 , where Raman spectra show the largest 
modifications after the complexation process, we con- 
clude that the geometry of the IMC-/3CD complex is the 
one shown in Fig. EI characterized by the 4-chlorobcnzoyl 
unit inserted into the cavity of /3CD through its larger 
rim, as early proposed by Fronza et al.pi^ We have also 
computed the density of vibrational states of free IMC, of 
/3CD, and of their inclusion complex, and compared it to 
the corresponding experimental Raman spectra. Though 
it is difficult to reach complete quantitative agreement 
between experiment and simulations, the latter explain 
the shift to lower frequencies of the most intense Ra- 
man peak, whose eigenvector is localized on the atoms of 
the amide C=0 group; this shift is due only to the inclu- 
sion process, whereas dimerization processes or structural 
changes in the solid state seem not to be involved. This 
multi-technique, experimental and numerical approach, 
appears very promising for the study of other inclusion 
complexes as well. 



VI. ACKNOWLEDGMENTS 

We would like to thank G. Mariotto for making the 
micro- Raman apparatus available, A. Sterni for technical 
assistance in the MS analysis and F. Pederiva for useful 
discussions. 



1 S. Santesson, J. Johansson, L.S. Taylor, L. Levander, S. (2003). 
Fox, M. Sepaniak and S. Nilsson, Anal. Chem. 75, 2177 



8 



2 G. Fini, J. Raman Spectrosc. 35, 335 (2004). 

3 S. Li and W.C. Purdy, Chem. Rev. 92, 1457 (1992). 

4 A.M. Amado, A.M. Moreira da Silva, P.J.A. Ribeiro-Claro 
and J.J.C. Teixeira-Dias, J. Raman Spectrosc. 25, 599 
(1994). 

5 P.J.A. Ribeiro-Claro, A.M. Amado, J.J.C. Teixeira-Dias, 
J. Raman Spectrosc. 27, 155 (1996). 

6 R. Spivey, R.L. Swofford, Appl. Spectrosc. 53, 435 (1999). 

7 T. Iliescu, M. Baia and V. Miclaus, Euro J. Pharm. Sci. 
22, 487 (2004). 

8 S.H. Choi, S.Y. Kim, J.J. Ryoo, J.Y. Park and K.P. Lee, 
Anal. Sci. 17, 1785 (2001). 

9 E. Redenti, M. Zanol, P. Ventura, G. Fronza, A. Comotti, 
P. Taddei and A. Bertoluzza, Biospectrosc. 5, 243 (1999). 

10 G.S. Zhang, C.Q. Tu, G.Y. Zhang, G.B. Zhou and W.L. 
Zheng, Leuk. Res. 24, 385 (2000). 

11 V. Andronis and G. Zografi, J. Non-Cryst. Solids 271, 236 
(2000). 

12 A.M.C. Myles, D.J. Barlow, G. France and M.J. Lawrence, 
Biochim. Biophys. Acta 1199, 27 (1994). 

13 G. Fronza, A. Mele, E. Redenti and P. Ventura, 
J. Org. Chem. 61, 909 (1996). 

14 L.S. Taylor and G. Zografi, Pharm. Res. 14, 1691 (1997). 

15 M. Fleischman, P.J. Hendra and A.J. McQuillan, 
Chem. Phys. Lett 26, 123 (1974). 

16 D.L. Jeanmaire and R.P.V. Duyne, J. Electroanal. Chem. 
84, 1 (1977). 

17 M.G. Albrecht and J.A. Creighton, J. Am. Chem. Soc. 99, 
5215 (1977). 

18 M.J. Frisch et al., Gaussian 03 (Revision A. 10), Gaussian, 
Inc., Pittsburgh, PA (2003). 

19 C. Lee, W. Yang, and R.G. Parr, Phys. Rev. B 37, 785 
(1988). 

20 A.D. Becke, J. Chem. Phys. 98, 5648 (1993). 



21 P.C. Hariharan and J.A. Pople, Theor. Chim. Acta 28, 213 
(1973). 

22 M.W. Wong, Chem. Phys. Lett. 256, 391 (1996). 

23 J.A. Pople, A.P. Scott, M.W. Wong and L. Radom, 
Isr. J Chem. 33, 345 (1993). 

24 DA. Pearlman, D.A. Case, J.W. Caldwell, W.S. Ross, T.E. 
Cheatham, III, S. DeBolt, D. Ferguson. G. Seibel and P. 
Kollman, Comp. Phys. Commun. 91, 1 (1995). 

25 J. Wang, R.M. Wolf, J.W. Caldwell, P.A. Kollman, D.A. 
Case, J. Comput. Chem. 25, 1157 (2004). 

26 P. Vassilev, C. Hartnig, M.T.M. Koper, F. Frechard, and 
RA. van Santen, J. Chem. Phys. 115, 9815 (2001). 

27 H.M. Berman, J. Westbrook, Z. Feng, G. Gilliland, T.N. 
Bhat, H. Weissig, I.N. Shindyalov, P.E. Bourne, Nucleic 
Acids Res. 28, 235 (2000). 

28 AW. Schuettelkopf and D.M.F. van Aalten, Acta Crystal- 
logr. D60, 1355 (2004). 

29 A. Jakalian, D.B. Jack, and C.I. Bayly, J. Comput. Chem., 
23 1623 (2002). 

30 C.S. Lu, C.J.Hu, Y. Yu and Q.J. Meng, 
Chem. Pharm. Bull. 48, 56 (2000). 

31 A.S. Lee and Y.S. Li, Spectrochim Acta, Part A 52, 16173 
(1996). 

32 B. Rossi, P. Verrocchio, G. Viliani, G. Scarduelli, I. 
Mancini, G. Guella and F. Rossi, to be published. 

33 X.Chen, K.R. Morris, U.J. Griesser, S.R. Byrn and J.G. 
Stowell, J. Am. Chem. Soc. 124, 15012 (2002). 

34 A. Jubert, N.E. Massa, L.L. Tevez and N.B. Okulik, 
Vib. Spectrosc. 37, 161 (2005). 

35 H.J.C. Beredsen, J. P.M. Postma, W.F. van Gunsteren, A. 
Di Nola and J.R.Haak, J. Chem. Phys 81, 3684 (1984). 

36 A. Anderson, The Raman effect, Vol. 1, Marcel Dekker, 
New York (1971). 



